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We study the effect of isospin degree of freedom on the balance energy (E M ) as well as 
its mass dependence throughout the mass range 48-270 for two sets of isobaric systems 
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ones. Our findings reveal the dominance of Coulomb repulsion in isospin effects on E& a ; as 
well as its mass dependence throughout the range of the colliding geometry. Our results 
also indicate that the effect of symmetry energy and nucleon-nucleon cross section on 
is uniform throughout the mass range and throughout the colliding geometry. We also 
present the counter balancing of nucleon-nucleon collisions and mean field by reducing 
the Coulomb and the counter balancing of Coulomb and mean filed by removing the 
nucleon-nucleon collisions. 
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1 Introduction 



With the availability of radioactive ion beam (RIB) facilities at Cooler Storage Ring 
(CSR) (China) [1], the GSI Facility for Antiproton and Ion beam Research (FAIR) [2], 
RIB facility at Rikagaku Kenyusho (RIKEN) in Japan [3], GANIL in France [4], and 
the upcoming facility for RIB at Michigan State University [5] one has a possibility to 
study the properties of nuclear matter under the extreme conditions of isospin asymmetry. 
Heavy-ion collisions induced by the neutron rich matter provide a unique opportunity to 
explore the isospin dependence of in-medium nuclear interactions, since isospin degree of 
freedom plays an important role in heavy-ion collisions through both the nuclear matter 
equation of state (EOS) as well as via in-medium nucleon-nucleon (nn) cross section. 

After about three decades of intensive efforts in both nuclear experiments and theo- 
retical calculations, equation of state for isospin symmetric matter is now relatively well 
understood [6]. The effect of isospin degree of freedom on the collective transverse in-plane 
flow as well as on its disappearance [7] (there exists a particular incident energy called 
balance energy (E ba i) or energy of vanishing flow (EVF) at which transverse in-plane flow 
disappears) has been reported in the literature [8-10], where it was found that neutron-rich 
systems have higher E ba i compared to neutron-deficient systems at all colliding geometries 
varying from central to peripheral ones. The effect of isospin degree of freedom on E bal was 
found to be much more pronounced at peripheral colliding geometries compared to central 
ones. As reported in the literature, the isospin dependence of collective flow as well as its 
disappearance has been explained as a competition among various reaction mechanisms, 
such as nucleon-nucleon collisions, symmetry energy, surface property of the colliding nu- 
clei, and Coulomb force. The relative importance among these mechanisms is not yet clear 
[8]. In recent study, we [11] confronted theoretical calculations (using isospin-dependent 
quantum molecular dynamics (IQMD) model [12]) with the data at all colliding geome- 
tries and were able to reproduce the data within 5% on average at all colliding geometries. 
Motivated by the good agreement of the calculations with data, two of us [13] (in order to 
explore the relative importance among above mentioned reaction mechanisms in isospin 
effects on the E ba i) calculated the E ba i throughout the mass range for two sets of isotopic 
systems with N/Z =1.16 and 1.33. The isotopic pairs were chosen so that the effect of 
the Coulomb repulsion is the same for a given pair. The choice of the above N/Z was 
taken because the percentage difference A N/Z(%) = ^-^\m^r^ — xlOO between the 
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given pair is about 15% which is same as in Ref. [9, 11]. Based on the results of Eb a / for 
isotopic pairs we concluded in Ref. [13] that mass dependence effects seem to dominate 
the isospin effects (consisting of isospin-dependent cross section, symmetry energy, surface 
properties). In our recent results [14] we find that collective flow for isotopic pairs with 
large difference between N/Z is sensitive to the symmetry energy. Coulomb repulsion will 
be same for isotopic pair throughout the mass range. The comparison of E& a / for isotopic 
pairs gave us the hint that the Coulomb repulsion could be dominant in isospin effects 
in collective flow as well as its disappearance [8,9, 11]. Therefore Gautam and Sood [13] 
studied the isospin effects on the throughout the mass range 48-350 for two sets of 
isobaric systems with N/Z = 1.0 and 1.4 at semi central colliding geometry. These results 
showed that the difference between the Eu for two isobaric systems is mainly due to the 
Coulomb repulsion. It was also shown that Coulomb repulsion dominates over symmetry 
energy. These findings also indicated towards the dominance of the Coulomb repulsion in 
larger magnitude of isospin effects in Ef, a i at peripheral collisions. Here we aim to extend 
the study over full range of colliding geometry varying from central to peripheral ones. 
Section 2 describes the model in brief. Section 3 explains the results and discussion and 
Sec. 4 summarizes the results. 

2 The model 

The present study is carried out within the framework of IQMD model [12]. The IQMD 
model treats different charge states of nucleons, deltas, and pions explicitly, as inherited 
from the Vlasov-Uehling-Uhlenbeck (VUU) model. The IQMD model has been used suc- 
cessfully for the analysis of a large number of observables from low to relativistic energies. 
The isospin degree of freedom enters into the calculations via symmetry potential, cross 
sections, and Coulomb interaction. 

In this model, baryons are represented by Gaussian-shaped density distributions 

1 1 2L 

fifrM = ^2 ex p(-[ f - f ^ 2 ^ x e M-W-Pi(t)} 2 -^) (i) 

Nucleons are initialized in a sphere with radius R = 1.12 A 1 / 3 fm, in accordance with 
liquid-drop model. Each nucleon occupies a volume of h 3 , so that phase space is uni- 
formly filled. The initial momenta are randomly chosen between and Fermi momentum 
(Pf)- The nucleons of the target and projectile interact by two- and three-body Skyrme 



3 



forces, Yukawa potential, Coulomb interactions, and momentum-dependent interactions. 
In addition to the use of explicit charge states of all baryons and mesons, a symme- 
try potential between protons and neutrons corresponding to the Bethe-Weizsacker mass 
formula has been included. The hadrons propagate using Hamilton equations of motion: 

dfl d(H}_ dpi d(H) 



dt dpi ' dt d?i 



(2) 



with 



(H) = (T) + (V) 

i 1 i j>i 

xfj(r',p',t)drdr' dp dp'. (3) 



The baryon potential V j , in the above relation, reads as 



v^(f ' - f) = vii + V2 , + V l J , + V ij ,- + V ij 

V ) Skyrme 1 r Yukawa 1 y Coul 1 y mdi 1 y sym 

= M(f -f)+t 2 5{f' -f)p'- l C-^)} 
^ exp(l(r'-f)l/jx) , Z^e 2 

+t 3 — 7TT=T, -^i / \ r 



(|(f'-f)|//i) \(f'-f)\ 

+t 4 \n 2 [t 5 (p' -p) 2 + l]5(f ' -r) 

+t 6 -T 3i T 3j S(r i '-r j ). (4) 
Qo 

Here Zi and Zj denote the charges of ith and jth baryon, and T 3i and T 3 j are their 
respective T 3 components (i.e., 1/2 for protons and —1/2 for neutrons). The parameters fi 
and ti,....,tQ are adjusted to the real part of the nucleonic optical potential. For the density 
dependence of the nucleon optical potential, standard Skyrme type parametrization is 
employed. The momentum-dependence V^ di of the nn interactions, which may optionally 
be used in IQMD, is fitted to the experimental data in the real part of the nucleon optical 
potential. It is worth mentioning that the Gaussian width which describes the interaction 
range of the particle depends on the mass of the system in IQMD, since each nucleus 
shows maximum stability for a particular width as shown in Ref. [11, 12] . For eg, width 
of Ca is 4.16 fm 2 and for Au is 8.33 fm 2 and Gaussian width varies between this mass 
range. 
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3 Results and discussion 



For the present study, we simulate several thousands events of each reaction at incident 
energies around Eu in small steps of 10 MeV/nucleon. In particular, we simulate the reac- 
tions 24 Mg+ 24 Mg, 58 Cu+ 58 Cu, 72 Kr+ 72 Kr, 96 Cd+ 96 Cd, 120 Nd+ 120 Nd, 135 Ho+ 135 Ho, hav- 
ing N/Z = 1.0 and reactions 24 Ne+ 24 Ne, 58 Cr+ 58 Cr, 72 Zn+ 72 Zn, 96 Zr+ 96 Zr, 120 Sn+ 120 Sn, 
and 135 Ba+ 135 Ba, having N/Z = 1.4, respectively, in the whole range of colliding geometry. 
The colliding geometry is divided into four impact parameter bins of 0.15 < b < 0.25 (BIN 
1), 0.35 < b <0.45 (BIN 2), 0.55 < b < 0.65 (BIN 3), and 0.75 < b < 0.85 (BIN 4), where 
b = h/h max . Here N/Z is changed by keeping the mass fixed. We use anisotropic standard 
isospin- and energy-dependent nn cross section a = 0.8 cf NN ^ ree . The details about the 
elastic and inelastic cross sections for proton-proton and proton-neutron collisions can be 
found in Ref. [12]. The cross sections for neutron-neutron collisions are assumed to be 
equal to the proton-proton cross sections. 

We also use soft equation of state along with momentum-dependent intercations 
(MDI). The results with the above choice of equation of state and cross section were 
in good agreement with the data [11, 13]. The reactions are followed until the transverse 
flow saturates. The saturation time varies form 100 fm/c for lighter masses to 300 fm/c for 
heavier masses. For transverse flow, we use the quantity " directed transverse momentum 
(pf r y which is defined as [15, 16] 

1 A 

(^^E^i/WKW, (5) 

i=l 

where y{i) is the rapidity and p x (i) is the momentum of i th particle. The rapidity is 
defined as 

where E(i) and p z (i) are, respectively, the energy and longitudinal momentum of i th 
particle. In this definition, all the rapidity bins are taken into account. A straight line 
interpolation is used to calculate E 6ai . 

Figure 1 displays the mass dependence of for four impact parameter bins. The 
solid (open) green circles indicate E^ for systems with lower (higher) N/Z. Lines are 
power law fit oc A r . Left (right) panels are for mass dependence of E^ when we include 
(exclude) A = 48. First we discuss the left panels for all the four bins. E bat follows a 
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Table 1: The values of T\ and t\a for BIN 1 to BIN 4 for calculations with Coulomb 
potential with and without A=48. 



b/b ma2: 


Tl 


Tl.4 


With A=48 


Without A=48 


With A=48 


Without A=48 


BIN 1 


-0.36± 0.01 


-0.38± 0.03 


-0.33± 0.01 


-0.31± 0.02 


BIN 2 


-0.50± 0.01 


-0.54± 0.03 


-0.45± 0.01 


-0.48± 0.01 


BIN 3 


-0.70± 0.03 


-0.83± 0.07 


-0.56± 0.02 


-0.57± 0.07 


BIN 4 


-0.93± 0.08 


-1.14± 0.25 


-0.70± 0.01 


-0.70± 0.03 



Table 2: Same as table 1 but for calculations with Coulomb potential reduced by a factor 
of 100. 



b/b max 


Tl 


Tl.4 


With A=48 


Without A=48 


With A=48 


Without A=48 


BIN 1 


-0.25± 0.02 


-0.17± 0.02 


-0.26± 0.03 


-0.14± 0.01 


BIN 2 


-0.25± 0.02 


-0.17± 0.02 


-0.28± 0.02 


-0.19± 0.02 


BIN 3 


-0.33± 0.02 


-0.24± 0.03 


-0.33± 0.02 


-0.26± 0.01 


BIN 4 


-0.41± 0.02 


-0.31± 0.03 


-0.41± 0.02 


-0.29± 0.02 



power law behavior oc A r for both N/Z = 1 and 1.4 (r being labeled as Ti.o and t xa for 
systems having N/Z = 1 and 1.4, respectively) at all colliding geometries. There are small 
deviations from power law behavior for heavy mass systems with N/Z = 1 at peripheral 
colliding geometry. Isospin effects are clearly visible for all the four bins as neutron-rich 
system has higher E^ai throughout the mass range in agreement with the previous studies 
[8,9,11]. 

The magnitude of isospin effects increases with increase in mass of the system at all 
colliding geometries. The effect is much more pronounced at larger colliding geometries. 
One can see that the difference between ri. and t\a increases with increase in the impact 
parameter (green circles). In Ref. [13], Gautam and Sood studied the isospin effects on 
the mass dependence of E 6ai for BIN 2. There they reduced the Coulomb potential by a 
factor of 100 and showed that the Coulomb repulsion plays dominant role over symmetry 
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System size (A) 



Figure 1: (Color online) Left (right) panels displays E M as a function of combined mass of 
system for different impact parameter bins with (without) A = 48. Solid (open) symbols 
are for systems having N/Z = 1.0 (1.4). Circles (diamonds) are for calculations with 
full (reduced) Coulomb. Lines are power law fit cc A T . r values without errors for full 
Coulomb calculations are displayed in figure. The detailed values of r are given in the 
table 1 and 2. 
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Figure 2: (Color online) Upper (lower) panels display r as a function of reduced impact 
parameter for full (reduced) Coulomb. Values are plotted at upper limit of each bin. 
Symbols are explained in text. Lines are only to guide the eye. 



energy in isospin effects on as well as its mass dependence at semi central colliding 
geometry (BIN 2). Since here we plan to extend that study over a full range of colliding 
geometry, so here also we reduce the Coulomb potential by a factor of 100 and calculate 
the E bal throughout the mass range at all colliding geometries. Solid (open) diamonds 
represent E& a j calculated with reduced Coulomb for systems with lower (higher) neutron 
content. Lines are power law fit oc A r . Interestingly, we find that: 

(a) the magnitude of isospin effects (difference in E& a j for a given pair) is now nearly same 
throughout the mass range which indicates that the effect of symmetry energy is uniform 
throughout the mass range. This is true for all the colliding geometries. This is supported 
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0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

b/b 

max 

Figure 3: (Color online) as a function of impact parameter for different system masses. 
Various symbols have same meaning as in Fig.l. Lines are only to guide the eye. 

by Ref. [17] where Sood and Puri studied the average density as a function of mass of the 
system (throughout the periodic table) at incident energies equal to E ba[ for each given 
mass. There they found that although both E^ai and average density follows a power law 
behavior oc A T , decreases more sharply with the combined mass of the system (with 
r = -0.42), whereas the average density (calculated at incident energy equal to E ba ;) is 
almost independent of the mass of the system with r = -0.05. It is worth mentioning here 
that the trend will be different at fixed incident energy in which case density increases 
with increase in the mass of the system [18, 19]. We also note that the power law with 
reduced Coulomb is now absolute for N/Z = 1 (solid diamonds fig. 1 (d)). This indicates 
that the deviations from power law behavior for heavy mass systems (with N/Z = 1) are 
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due to the dominance of Coulomb repulsion. 

(b) one can also see that the enhancement in E^ (by reducing Coulomb) is more in 
heavier systems as compared to lighter systems for all colliding geometries. The effect is 
more pronounced at higher colliding geometries. 

(c) throughout the mass range at all colliding geometries, the neutron-rich systems have 
a decreased E& a ; as compared to neutron-deficient systems when we reduce the Coulomb. 
This trend is quite the opposite to the one which we have when we have full Coulomb. 
This (as explained in Ref. [13] also) is due to the fact that the reduced Coulomb repulsion 
leads to higher E^. As a result, the density achieved during the course of the reaction 
will be more due to which the impact of the repulsive symmetry energy will be more in 
neutron-rich systems, which in turn leads to a decreased for neutron-rich systems. 

Now to discuss how the value of r changes if we exclude lighter systems (right panels 
in fig. 1), we display in fig. 2 the variation of r as a function of impact parameter. Solid 
(open) symbols are for full (reduced) Coulomb. Triangles (pentagons) represent r with 
(without) A = 48. For full Coulomb (upper left panel), T\ A increases with increase in 
impact parameter but the increase is independent of inclusion/exclusion of lighter mass. 
On the other hand, T\ increases drastically with impact parameter (upper right panel). 
Moreover, the increase in T\ (with impact parameter) is more sharp when we take into 
account only heavier masses (pentagons) as compared to when lighter mass systems are 
also included (triangles). Thus as we have discussed in Ref. [13] also that if we take 
into account only heavier systems, the value of r is more as compared to if we include 
lighter systems as well. Here we see enhancement of this effect with increase in colliding 
geometry showing the increased role of Coulomb repulsion in mass dependence of E ba i at 
high impact parameter. In lower left panel, for reduced Coulomb, we see the values of T\ 
and T1.4 increases much less sharply with increase in colliding geometry and this increase is 
almost independent of N/Z of the system which shows that the effect of isosp in-dependent 
cross section and symmetry energy on the mass dependence of E M is independent of N/Z 
throughout the colliding geometry. The same is true if we consider only heavier systems 
as well (lower right panel). As we have seen in fig. 1, the effect of symmetry energy is 
uniform throughout the mass range at all colliding geometries, this means that the effect 
of cross section must also be nearly same through the mass range and colliding geometry. 
We will come to this point later. It is worth mentioning that since IQMD reproduces the 
data nicely [11] so the experimental values of tau are not expected to differ much from 
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Figure 4: (Color online) at central and peripheral colliiding geometries for A = 116 
(upper panel) and A = 240 (lower panel) with no nucleon-nucleon collisions (Hexagons). 
Circles represent the calculations with collisions. 



the values of tau given in table 1. 

In fig. 3a, 3b, and 3c, we display E« as a function of b for masses 116, 192, and 240, 
respectively, for both full and reduced Coulomb. Symbols have the same meaning as in fig. 
1. For full Coulomb (green circles), for all the masses at all colliding geometries, system 
with higher N/Z has larger Em in agreement with previous studies [8,9,20]. Moreover, 
the difference between for a given mass pair, increases with increase in colliding 
geometry. This is more clearly visible in heavier masses. Also for N/Z = 1.4, E& a ; increases 
with increase in impact parameter in agreement with Ref. [20]. This is due to the 
decreased participant zone in peripheral collisions which decreases the amount of repulsive 
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Figure 5: (Color online) The time evolution of collision rate ( ?? u ) for various system 
masses at 50 MeV/nucleon for BIN 1 and BIN 4. Various lines are explained in the text . 

nn collisions and therefore higher energy is required to overcome the attractive nuclear 
force. This effect is less pronounced in heavier systems since even at peripheral geometry 
significant number of nn collisions will occur. Moreover, the effect of stronger Coulomb 
repulsion in heavier systems will increase with colliding geometry since it will push more 
number of nucleons in the transverse direction away from participant zone. However for 
N/Z = 1, increase in with impact parameter is true only for lighter mass system such 
as A = 116. For heavier masses E^ infact begins to decrease with increase in impact 
parameter in contrast to the previous studies [8,9, 11,20]. However, when we reduce the 
Coulomb (by a factor of 100 (diamonds)), we find that: 
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(i) Neutron-rich systems have a decreased Eu as compared to neutron-deficient sys- 
tems as mentioned previously also. This is true at all the colliding geometries throughout 
the mass range. This clearly shows the dominance of Coulomb repulsion over symmetry 
energy in isospin effects throughout the mass range at all colliding geometries. 

(ii) The difference between E bai for systems with different N/Z remains almost constant 
as a function of colliding geometry which indicates that the effect of symmetry energy is 
uniform throughout the range of b as well. This also shows that the large differences in 
Eb a i values for a given isobaric pair are due to the Coulomb repulsions. 

(iii) In heavier systems, at high colliding geometry, the increase in E ba i is more in 
systems with N/Z =1 compared to N/Z = 1.4 which shows the much dominant role of 
Coulomb repulsion at high colliding geometry. 

To see the relative contribution of Coulomb repulsion and cross section in lighter and 
heavier systems, we switch off the collision term and calculate the for A = 116 and 
240 at two extreme bins. The results are displayed in fig. 4. Hexagons represent the 
calculations without collisions. Other symbols have same meaning as in fig. 1. We find 
that at a given impact parameter Eu increases by large magnitude for both systems which 
shows the importance of collisions. The increase is of the same order for both the masses 
at both impact parameter bins indicating the same role of cross section for both lighter 
and heavier masses as we have expected in discussion of fig. 2. This is supported by 
Ref. [17, 18], which shows that since mean field is independent of the mass of the system, 
so one needs the same amount of collisions to counter balance the mean field in both 
lighter and heavier masses. Moreover, the same order of increase of E ba i at central and 
peripheral colliding geometry (when we switch off the collisions) indicates the importance 
of collisions at high colliding geometry as well. The effect that E ba i decreases with increase 
in impact parameter (due to dominance of Coulomb) for heavy mass systems with N/Z 
= 1 (fig. 3 lower panel) does not appear here for lighter and heavier masses. Therefore 
in fig. 3, the reduced Coulomb allows one to examine the balance of nn collisions and 
mean field while in fig. 4 the removal of nn collisions allows one to examine the balance 
of Coulomb and mean field. 

As we have seen in fig. 4 that E ba[ is much higher than the actual E ba[ when we switch 
off the cross section, so to explore whether Coulomb repulsion affects the collisions in the 
E bat domain, we display in fig. 5 collision rate dN ^° u for A = 48 (upper panel), 116 (middle 
panel), and 240 (bottom panel) at incident energy of 50 MeV/nucleon. Solid (dashed) 
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lines represent Coulomb full (reduced) calculations. Higher (lower) peaks represent results 
for central (peripheral) impact parameter. We find that Coulomb decreases the collision 
rate in medium and heavier mass systems for central bins whereas for peripheral bins, the 
effect of Coulomb on collisions is only for heavier masses (lower panel). Comparing top 
and bottom panel we see that there are still significant number of collisions at peripheral 
geometry. 

4 Summary 

We have studied the isospin effects in the disappearance of flow as well as its mass de- 
pendence throughout the mass range 48-270 for two sets of isobaric systems with N/Z = 
1 and 1.4 in the whole range of colliding geometry. Our results clearly demonstrate the 
dominance of Coulomb repulsion in isospin effects on E ba[ as well as its mass dependence 
throughout the range of colliding geometry. The above study also shows that the effect 
of symmetry energy on and cross section is uniform throughout the mass range and 
colliding geometry. We have also presented the counter balancing of nn collisions and 
mean field by reducing the Coulomb and the counter balancing of Coulomb and mean 
filed by removing the nn collisions. 

This work has been supported by a grant from Indo-French Centre For The Promotion 
Of Advanced Research (IFCPAR) under project no. 4104-1. 
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